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ABSTRACT 


The goal of this research was to experimentally депоп- 
strate the feasibility of a cylindrical shock wave to serve 
as the major pressure recovery mechanism for either 
cylindrical gas dynamic lasers or chemical lasers. An 
analysis was conducted to determine the minimum distance 
between the nozzle exit cylinder and the shock wave location, 
Et $ required for power extraction. A cylindrical nozzle 
array was designed using the axisymmetric method of charac- 
teristics for a Mach number of 4,5 with a flexibility of 
other Mach numbers. A chamber incorporating end wall boundary 
layer bleed and ejectors for back pressure control was 
designed and built. The apparatus was run up to a presSure 
of 18 atmospheres. Three flow regimes were identified. A 
value of (а 06) / (55750) of 0.3 was attained, where г 15 
the shock wave location. Flow in the end wall nozzles 
consisted of a Coanda jet with a subsonic flow region 
adjacent to the wall; flow was not altered by boundary layer 


bleed. Except for the end wall nozzles a cylindrical shock 


wave was established. 
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I. INTRODUCTION 


At present there are three types of lasers that employ 
supersonic flow: chemical, CO electrical, and gas dynanic. 
All three are capable of producing large amounts of lasing 
power, but are restricted to use in platforms that are not 
limited by weight and volume. Due primarily to the large 
diffusers necessary to restore the flow to ambient pressure, 
large volume and weight are required to handle a given flow. 

A cylindrical geometry gas dynamic laser has advantages 
not found in the other supersonic lasers mentioned: 


(1) diffuserless operation 
(2) elimination of nozzle bending stresses 


(3) elimination of boundary layer corner interactions. 
Such a geometry would result in a compact system achieving 
lasing by rapid expansion in the radial direction rather 
than linearly as is presently done. The supersonic radial 
flow is then slowed to subsonic and raised to near ambient 
pressures by traveling through a cylindrical shock wave 
sufficiently downstream of the nozzles to allow for a 
reasonably large lasing cavity. A cylindrical shock wave 
eliminates the need for diffusers, thereby reducing the 
weight and volume of the device considerably. 

Supersonic conditions are achieved from flow traveling 
through axially symmetric nozzle es iius eliminating 
nozzle bending stresses found in linear gas dynamic lasers. 


A schematic of a cylindrical GDL is shown in figure (1). 
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Figure (2) shows schematic diagrams of two types of 
linear GDL's. Boundary layer buildup on the four walls of 
the lasing cavity interact at the corners, creating increased 
possibilities of flow separation and non-homogeneous conditions 
in the lasing medium. By having a flat plate at each end of 
the nozzle stack, the cylindrical geometry would eliminate 
such interaction. Since the radial flow experiences 
increasing area as it moves outward, the boundary layer buildup 
on each plate is less than in the linear case, 

A cold flow test of a number of nozzle rings was conducted 
to establish a cylindrical shock wave and to determine whether 
such an apparatus is feasible for use in a cylindrical gas 
dynamic laser, 

A theoretical power output example is included in section 
II of this report. Section III contains the development of 
the nozzle feed, the nozzle ring geometry and the ejector and 
boundary layer bleed systems. Experimental results are 


described in section IV, 
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II. CYLINDRICAL LASER GEOMETRY 


A. PREVIOUS GDL'S 

Linear gas dynamic lasers are currently in use and have 
successfully produced laser beams in the order of a megawatt 
of power. These devices are quite large and heavy, thus 
limiting their use in operating environments (such as an 
aircraft) where such parameters are limited. These limita- 
tions are due primarily to the large diffusers necessary to 
return the supersonic flow to ambient pressure. 

A cylindrical GDL could reduce the weight and volume 
necessary to handle a given flow by expansion to supersonic 
velocities in the radial direction and returning to ambient 
pressure by means of a cylindrical shock wave. Other advan- 
tages include elimination of bending stresses and reduced 
boundary layer interferences. 

Ortwerth [Ref. (1)] designed a cylindrical source flow 
device, but retained diffusers similar to those used in 
linear GDL's. The diffusers were aligned radially as shovn 
in figure (3). Such a design reduces the volume required of 
such a device by having an annular lasing cavity rather than 
a linear one. Retaining the diffusers does not allow for 
the substantial reductions in weight and volume that a 
cylindrical shock wave would facilitate. 

Sponsored by the Air Force Weapons Laboratory, the AVCO 


Everett Research Laboratorv [Ref. (2)] conducted an extensive 
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theoretical program to establish the feasibility of 

cylindrical source flow for a cylindrical GDL. Included 

in the program was a cold flow test of one nozzle ring which 

established supersonic flow at the predicted mach number but 

failed to create a cylindrical shock wave. Results showed. 

that homogeneity of the lasing medium can be made 

acceptable. Separation of the boundary layer interfered with 

the establishment of a cylindrical shock wave, and it was 

recommended that some type of boundary layer control be used. 
Figure (4) shows the assembly used in the AVCO Everett 


Fest, 


В. THEORETICAL POWER OUTPUT 

In a linear GDL the diffusers can be placed far enough 
downstream so that a large lasing cavity is assured. In 
a cylindrical GDL the width of the lasing cavity is deter- 
mined by the shock location. To extract all the power 
available from a given flow, the cylindrical shock must not 
Occur too near the nozzle exit. Figure (5) shows the nozzle 
stack with the lasing mirrors. 

Assuming a Gaussian intensity distribution on each mirror, 
the intensity I is given by 
-r?c (1) 


I(r) = leak? 
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Satisfying the boundary conditions and re-defining r yields 


а е e kilowatts (2) 


P cm 
where r, and r are shown in figure (A-1) in Appendix A. 


Integrating over the mirror area gives the total power 


Tm 


Ср 12 ar (3) 
٣ے‎ E Y 


An alternate method of finding the power available from 


the flow is 
Power = (mass flow rate) x (specific energy) 
er 
E 2 те 10.80 (0/0) (а/а) H. x (E) kilowatts (4) 


where subscripts "e" and "o" refer to exit and stagnation 


conditions respectively.  Equating equations (3) and (4) gives 


E 
m 


r-r 2 _ 
In ы Y pl V2 1° ак = Ге 202240 (p/p) (a,/a,) M. E 
e e (5) 


where the only unknown is the outer radius of the mirror, 


E 
m 


Zu 





Performing the integral in equation (5) and solving for 


г тсүг + с | (6) 


ar, Іра (а/а) М. Е (р/р) 


72 І, erf (V2) 


Using the values for the constants listed in Appendix A, 


equation (6) reduces to 


г #У2,13 + .2692 L (7) 


m 


where L is the overall length of the nozzle stack. Figure 
(6) shows the minimum shock location for increasing nozzle 
stack length. 

For six nozzle rings, the outer radius of the mirrors 
iS7 inches. Using this value in equation (3), the 
power is found to be 32 kilowatts. Details of the theoretical 


power available calculations are shown in Appendix A, 
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111. EXPERIMENTAL DEVELOPMENT 


Ge DESCRIPTION OF EXPERIMENT 

Reference 2 is a final report of analytical and experi- 
mental programs conducted by the AVCO Everett Research 
Laboratory to determine the feasibility of using a cylindrical 
geometry for a gas dynamic laser. The Naval Weapons Center, 
China Lake, California, expressed an interest in the concept 
of a cylindrical GDL and requested that research be conducted 
to determine whether a cylindrical shock wave could be 
established by incorporating the recommendations of the AVCO 
Everett report. Funding of this project was provided by 
Naval Air Systems Command, Code 320, 

Before determining the actual shape of the nozzle rings 
(see figures (8) ana (9)), it was first necessary to calculate 
the throat height, that is, the spacing between each nozzle 
ring in the stack. This spacing is determined by the mass 
flow that could be supplied to the nozzle stack, which is 
limited by the feed holes. See figure (14). 

Once the throat height was determined, the exit to throat 
area ratio, and thus the exit Mach number, was found from the 
geometry of the nozzle rings. The axisymmetric method of 
characteristics found in Ref. 3 was then used to determine 
the contour of the nozzle walls. Once the equation describing 
the wall shape was calculated, the design was sent to the 


machine facility at NWC, China Lake, where the nozzle rings 
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were fabricated. Figures (8), (9), (10), (12) and (14) were 
provided by NWC, China Lake, 

To prevent separation of the boundary layer at the shock 
location, suction was used to bleed off the low energy flow 
in the lower portion of the boundary layer. Extensive 
analytical investigation of supersonic boundary layer bleed 
was conducted and is found in section III-D of this report. 

The loss in total pressure across a normal shock of the 
strength anticipated resulted in total pressure downstream 
of the shock (using a reservoir pressure of ten atmospheres) 
to be less than one atmosphere,  Ejector nozzles were 
necessary to lower the back pressure seen by the flow exiting 
the nozzle rings and thus allow them to start. Design of 
the ejector apparatus is found in section III-E of this 
Report. 

Figure (7) shows the assembled experimental apparatus 
with dimensions. The design allowed for testing of up to 
eight nozzles at the design Mach number, The assembly shown 
was attached to the free jet located in building 230 of the 
Naval Postgraduate School. Air was supplied to the free jet 
from a 2400 cubic foot tank capable of pressure up to three 
hundred pounds per square inch. The moisture was removed 
from the air before it entered the tank by means of a 
hydryer. A model 06231 Lammert Industries vacuum pump was 
used to evacuate a 200 cubic foot tank to provide suction 


for boundary layer bleed. 
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DISASSEMBLED NOZZLE STACK (COURTESY NWC, CHINA LAKE) 


FIGURE (9) 





В. NOZZLE FEED SYSTEM 
1. Throat Height 

It was found, as expected, that the throat height 
was the most critical dimension of the entire design. A 
tively large throat height is desirable in that the 
ratio of the nozzle exit area to the throat area is kept 
small, thus keeping the exit Mach number reasonably low. 
This would have the advantages of ease in starting and 
smaller total pressure loss across the shock. 

A large throat height also requires large mass flow 
requirements, and the amount of air that can be fed to the 
nozzle rings is limited by the inlet holes in the center of 
each ring. Figures (10) through (14) show the dimensions 
and photographs of the nozzle rings. The area of the feed 
holes must be greater than the combined throat areas of the 
nozzle ring stack to prevent choking of the flow at the feed 
holes. Thus if many nozzle rings are to be tested, an 
extremely small throat height must be used. This would 
result in a number of complications: 


(i) An excessively large A causing 


exit/*throat’ 
high supersonic or hypersonic flow at the nozzle 
exit and thereby creating a problem in starting 


and a large loss in total pressure across the shock. 


(11) Extreme heat concentrations at the throat (not a 
problem in cold flow testing). 

(iii) Difficulty in fabricating spacing washers of uniform 
thickness, causing unacceptable variations in the 


flow due to changes in the throat height. 
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FIGURE (13) DIMENSIONS OF NOZZLE RING (END VIEW) 
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END VIEW OF NOZZLE STACK 


FIGURE (14) 





The throat height must satisfy the continuity 


equation: 


where n is the number of nozzle rings. This can also be 
expressed as 
орви. А. y? = 6p*U*A* (1) 
In in ain : 
Rewriting equation (1) uSing the equation of state and 


solving for the throat height: 


Е J rep kiwa 
В = ет Мазар) (т. 
іп 
Using the value of inlet area previously obtained and using 
the isentropic relations for pressure and temperature in 
equation (2), the throat height in inches can be expressed 


as 


E 1092 и | (3) 


295 
(1+.2M,,) 


Tote factor of two in equation (1) results from the 
nozzle stack being fed from both ends. 
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Equation (3) shows that the throat height is a function of 
inlet Mach number only. Details of the throat height 
calculations are found in Appendix С. 

The exit Mach number is determined by the ratio of 
the exit area to the throat area, This ratio is found from 


the geometry of figure (11) to be 





Веста _ 27r, (htt) 
A 2nr,h 
ог 
А 


СЕВ. 5) 
ca Te n 
Also, assuming isentropic flow through the nozzle, the area 


ratio can be expressed as 


YL 
. gr (p м 21? р. (5) 
Equations (4) and (5) show that the exit Mach number is 
directly dependent upon the throat height. Since the throat 
height is a function only of the inlet Mach number, the exit 
Mach number is also a function of the inlet Mach number, 

The results of equations (3), (4) and (5) for various inlet 
Mach numbers are listed in Table (1). Figures (15) through 


(17) show the results of Table (1) in graphical form. 
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Table 1 Throat Height and Exit Mach Number for 


Various Inlet Mach Numbers. 


mier h (inches) МЕ 
„2 . 0179 5. б> 
„4 . 0334 4.95 
. 6 . 0446 4.62 
8 2051 4.47 
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A throat height of „05 inches was chosen to allow 
for the lowest exit Mach number that would not result in a 
choked inlet condition. From figure (17), a throat height 
of „05 inches results in an exit Mach number of 4,5. The 
axisymmetric method of characteristics was used to design 
the contour of the nozzle rings. These calculations are 
found in section III-C. 

Knowing the throat height between each nozzle ring, 


the mass flow rate can be determined. For each nozzle, 


= АЖТЖКАЖ 
Mhozzle en i 


Using the equation of state and rearranging yields 


DS 
Р (р) (49.1) 
O 


"n T* N 144 
O То 


ЕТО) 
— 
Substituing the appropriate values, the mass flow per 
nozzle is found to be 


Mm, = 1.08 lb m/sec . 


Total mass flow is 


Me otal 6 (my) 


т = 6.48 lb /Sec 
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2. Axial Variation of Mach Number and Static Pressure 
Arter the flow enters the nozzle stack, changes in 
Mach number and static pressure occur in the axial direction. 
Shapiro [Ref. 4] develops influence coefficients that relate 
changes in one flow condition to changes in another. 
Assuming constant molecular weight, the change in Mach 


number is related to the area change by 


Y-1 ,,2 
ам TA (6) 
М 1-М 


Appendix D shows the model used for the axial flow through 
the nozzle stack. The pressure change is also related to 


changes in area and is given by 


E A O (7) 
A 


Using the geometry of the nozzle stack, and rearranging 
equations (6) and (7), axial Mach number and pressure 


changes are found to be 


y-l 2 
Lt = My 
M; = М. - M; c 1 90950) (8) 
1-м. 
1 
апа 
y 9) 
ER Ви LZ (.06667)] . ( 
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A computer program was written using equations (8) and (9) 
to determine values of Mach number and pressure in the 
axial direction. Figures (18) and (19) are the results of 
this program in graphical form. A listing of the program 
and the results are found in Appendix D, along with details 
of the preceding development. 

3. Radial Variation of Pressure and Mach Number 

As the flow leaves the nozzles, it experiences 

increasing area as it moves out radially. This area 
increase results in increasing Mach number and a corresponding 
decrease in static pressure. 


From the nozzle ring geometry, 


me 2m st 
A/A* - urere 
ме). (10) 


Assuming isentropic flow: 


Y+1 
_ 1 2 ү-1 ,,2 2 (у-Т) 
А/А* = М F p = МЕЗІ (11) 
and 
КҮТЕ, t 
E ШЕ из 
A SS A) | (12) 
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Equations (10), (11) and (12) relate Mach number and 
pressure to radial position. Table (2) is a summary of 
pese three equations for different radii. Figures (20) 
and (21) are graphical representations of Mach number and 


pressure as functions of radius. 


С. NOZZLE RING DESIGN 
MX Nozzle Contour 
As was shown previously in Section III-üA, an exit 
Mach number of 4,5 was chosen for design. The curvature of 
the nozzle rings was determined by the use of the axisym- 
metric method of characteristics. 
актыла Turning 
In a two-dimensional nozzle the flow is turned 
is 


enrough (wv /2) degrees at the throat, where v 


final final 


the value of the Prandtl-Meyer function corresponding to the 
design exit Mach number. In a three-dimenSional case the 
flow expands from throat to tip (and therefore the Mach 
number increases) due solely to the increase in area, Thus 
to achieve the desired exit Mach number in the three dimen- 
sional case, the flow should be turned less than in the 
two-dimensional case, 

Assuming parallel plates, the ratio of the exit 
area to the throat area reduces to 
В - 


AS Benroat u 
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Table 2 Radial Variation of Mach Number and Pressure 


r (in.) n o stat PSR 
eo 4.50 SEL 
2.0 4.84 48.30 
2.5 5,089 30:07 
2,0 52239 302539 
ee) 5.479 PAS аи 
4,0 5. 674 18.85 
ше 5 5.869 14.92 
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Using a design Mach number corresponding to an area ratio 


of 


ll 


DUE CA) RE 
or 


A/A* 


tl 


О 


and then treating the initial turning as if it were a 
two-dimensional case should give approximate results (that 


is, an exit Mach number near 4,5). For 
A/A* = 10.4 , 
м ., = 3,96 
CATE 


which yields 


d 5523" 


exit 


Therefore the initial turning of the flow at the throat 
should be 
Vexit 


АӨ = —á——- 32.65? 


where 0 is measured from the flow direction entering the 


ehroat., 
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b. Axisymmetric Method of Characteristics 
Liepmann and Roshko [Ref. 3] develop expressions 
for finding the Prandtl-Meyer function (v) and the flow for 


direction (8) for an axially symmetric flow and are given by 


e nud d 
V3 = я (1 +52) + (901-6) 


1 | Sin 91 | sin 9, 
+ 5151п ц 2 4813 + Sin Wa —— Апоз! 
(13) 
and 
AS 
3 1572 250% 
1 | sin 6 | sin 9, 
IM у эп Ш м, 1 


(14) 


where all angles and subscripts are as shown in figure (22). 

Using equations (13) and (14) the nonsimple core 
of the characteristic mesh was calculated. Figure (23) shows 
the entire characteristic mesh. The initial turning at the 
sharp-edged throat was divided into eight segments as figure 
(23) illustrates. Section II of Appendix E contains tables 
of Mach number, Mach angle, flow angle and Prandtl-Meyer 
function for each point in the mesh. 

As the last of the right running characteristics 
from the throat reached the plane of symmetry (point 70), 


the Mach number was found to be 4,47, and the flow angle was 








|5 


Geo 


FIGURE (22) CHARACTERISTIC ELEMENT 
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90”, These flow conditions show that by reducing the 
design Mach number by the method previously described, and 
then treating the initial turning as a two-dimensional Васа 
accurate results will be obtained. 
с. Establishing Points on the Nozzle Wall 

The Tert running Characteristic from point 70 
is a straight line, and the flow conditions at any point on 
this characteristic will be the same as at point 70, that is, 
M= 4,47 and 8 = 90°, Points on the nozzle wall were deter- 
mined by calculating points along a right running character- 
istics in the -¿ direction starting from a point on the left 
running characteristic between points 70 and 76. 

Following a development similar to Ref, 2, 
equations (13) and (14) were modified to calculate points 
of the characteristic mesh proceeding in the -£ direction. 


The modified equations are 


E sde 
уҙ 5 %(У1%Уу2) + «(8, - 85) 


1 sin 0 Sin 8 
+ şl-sin Uu г. А61: + sin y, =. Апоз) 
(15) 
and 
"n "ES 1 
93 = xv, və) т VASE T ө.) 
1 | sin 9; sin Ө 
+ zl-sin u] cxi 2) Ur An) 3] 
(16) 
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where angles and subscripts are as shown in figure (E-1) 
found in Appendix E. Equations (15) and (16) were used to 
construct right running characteristics in the -E direction 
from points 72, 74, 76 and 78 (not shown) until the left 
running characteristics reflected from the centerline were 
met. Starting at the throat the wall was drawn at 32.65° 
(relative to the centerline) until point 19 was reached, 
where the wall angle was reduced to the flow angle that was 
calculated for point 19. This procedure was continued for 
points 32, 37 and 50. The last point (point 80) was the 
result of changing the angle of the wall to the estimated 
flow angle at that point. This line was extended until a 
radius of 1.5 inches was reached. Limiting the exit radius 
to this value was an effort to keep the exit Mach number 
as low aS posSible as was discussed previously. 

d. Wall Shape 

The r and z coordinates of the points where the 
wall slope changes, as well as the throat and the exit points, 
were fit with a fourth-degree polynomial using the Hewlett 
Packard model 9830-A computer and a polynomial regression 
program, The coordinates of the points and the coefficients 
of the fourth-degree polynomial are shown in Table (3). 
Figure (24) shows the nozzle contour from throat to exit. 
2. Trailing Edge 
As can be seen from figure (24), the nozzle contour 


is not tangent to the ring centerline at r = 1.5 inches. 
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Table 3 Description of Wall Contour 


A. Wall Point Coordinates from Method of Characteristics 


Point Number Соога, У Сока. 
(throat) 1.00 2025 
19 1.06 . 066 
52 Terl . 097 
37 1.19 „138 
50 12/29, 721413 
80 1,44 „215 
(exit) 2.50 226 


В. Coefficients of Fourth Degree Polynomial Fit to Points 





B(0) = -6,375 
B(1) = 18.144 
B(2) --19,773 
BI) ZS 9 
В (4) = -1.910 


R-Square = 0.99997 
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Rather than change the shape of the wall to a contour that 
would result in a sharp trailing edge, it was decided to 
maintain the design contour, This resulted in the flow 
exiting the nozzle at a 5° angle and a blunt trailing edge 
as shown in figure (25). 
3. Throat Shape 

A sharp-edged throat was used to obtain the highest 
degree of turning for the shortest length. Reference 5 
states that to be assured of uniform flow at the nozzle 
entrance, the radius of curvature of the inlet should be 
two to three times the throat height. Figure (26) shows the 
details of the nozzle inlet and throat. Due to the original 
geometry of the nozzle ring (before machining) the ratio of 
the radius of curvature to the throat height was limited to 
between 1.5 and 2,0. The shaded area in figure (26) 
represents the machining tolerance of the inlet contour. 
The center of curvature was displaced to the left of the 


sharp corner to assure horizontal tangency of the flow. 


р. EJECTOR NOZZLE DESIGN 

For a reservoir stagnation pressure of ten atmospheres 
and Mach-five flow, the total pressure downstream of the 
normal shock will be less than one atmosphere. To lower 
the back pressure and thus allow the cylindrical nozzles to 
start,ejector nozzles were designed and installed in the 
exhaust section of the apparatus as illustrated in figure 


(7). Using the continuity, momentum and energy equations, 
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FIGURE (26) DETAIL OF NOZZLE RING THROAT 
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solutions to the flow were found for different combinations 
of exhaust area and ejector nozzle diameter. A computer 
program was written to solve the three equations, and the 
results are tabulated in Appendix F. Figure (27) is a graph 
of these results. Any combination of nozzle diameter (Dl) 
and exhaust. width (D3) which results in a chamber pressure 
(P1) below the total pressure line will allow the cylindri- 
cal nozzles to start. A nozzle diameter of 1.75 inches and 
exhaust width of 6.5 inches were chosen. This combination 
results in a back pressure of „562 atmospheres. Complete 


details of this development are found in Appendix F. 


E. BOUNDARY LAYER BLEED 

Boundary layer bleed vas employed to keep the flow from 
separating at the shock location. By bleeding off the low 
energy air in the boundary layer the remaining flow will 
have sufficient energy to support the pressure rise associ- 
ated with a shock of the strength anticipated. A search of 
the literature relating to boundary layer failed to locate 
specific information on supersonic boundary layer bleed 
techniques, specifically, hole size and pressure ratio 
required to remove a certain fraction of the boundary layer. 
Since information could not be found, a model was developed 
to make such a prediction. 

1. Flow Model 

A model fluid: discussed by Stewartson [Ref. 7] was 

used in this development. The fluid is assumed to have а 


Prandtl number of unity and the viscosity is proportional to 
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E tatic temperature so that Пр = НР.” where subscript 1 
denotes the reference condition. Assuming the fluid has a 
Blasius profile and that stagnation temperature is constant 
through the boundary layer, the incompressible profile was 
transformed to a compressible profile by the method described 
in Ref. 7, 

For analytical ease, a Blasius profile was 


approximated by 


u/U = sin(An) 


where А = 54° and n is given by Ref. 7 as 


The reference length is denoted by X. The model shown in 
figure (28) was used to describe the flow. By applying a 
pesa plenum pressure (P,) lower than the freestream 
pressure (Pi), a Streamline (4) is assumed to stagnate on 
the downstream edge of a slot of width H. All flow below 
this streamline goes through the slot, which is assumed 
long enough to have parallel streamlines at the slot exit. 
By assuming isentropic flow in each streamtube below (4), 
it is possible to calculate flow conditions at the slot 
exit by specifying (Р./Р,). Using the continuity equation 


leads to a slot width, H. The complete set of equations 
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programmed for the Hewlett Packard model 9830 computer, 
and is listed in Appendix G. 
2. Results of Calculation 

Figure (29) illustrates the results of the computer 
program and plots the fraction of the boundary layer that 
enters the slot (Y/$) as a function of the pressure ratio 
(P, /P>) , with nondimensional slot width H/Š as a parameter. 
Table (4) gives the results for the flow conditions antici- 
pated for radial flow at a Mach number of 4,5, The symbols 
used are shown in figure (28) or have been given previously 
with the exception of во! which is the stagnation speed 
of sound at the edge of the boundary layer. 

The dot-dash-dot line in figure (29) represents the 
locus of maximum Y/ó for a given slot width and pressure 
ratio. 

3. Bleed Hole Design 
Using the values from table (4), (specifically H and 


X) the total bleed area was found to be 


Apleed = 2T (rg +X)H 


NSF) (00143) 


> e 
cc Ее ° 


Three hundred holes were drilled through each plate at a 


diameter of .141 inches as shown in figures (30) and (31). 
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FIGURE (29) PLOT OF FRACTION BOUNDARY LAYER BLEED, Y/6, 
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Table 4 
M = 4.5 
Y = 1.4 
U. = 2318 ft/sec 


6 = 0.00504 ft 


DE :00143 ft 


STREAMTUBE 
NUMBER n 
0 
1 
SE 
2 
32 
3 
48 
4 
. 64 
5 
. 80 


. 80 


41419 


1.58 


1595 


В. = 500,000 x = 0.25 ft 
A = 54° a = 1157 ft/sec 
ру = 200 psf pl 50 psf 
Ү/6 = 0.385 Н/6 = 0,283 
In ccu 1.82 E-4 slugs/sec 
AY AH M M. m —AH/N 
E-4 E-5 D 
£t ft slugs/sec 
.40 . 60 . 076 12506 0.74 0.15 
.40 шев “227 1,58 ee 0.47 
E39 3.01 ‚381 103 3.67 DOE 
. 38 3.98 2557 1.69 5.10 1.04 
597 4,75 .698 15 655) 1.28 
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Flow Conditions for Boundary Layer Bleed (from Ref. 
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FIGURE (30) BLEED HOLE DESIGN 
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BLEED HOLES WITH NOZZLE STACK 


FIGURE (31) 





The total area of the holes on one plate is larger than 


the calculated bleed area. This was done for three reasons: 


1) 


2) 


3) 


to bleed off a greater fraction of the boundary 
layer than calculated if needed (see figure (29)), 
to keep the same bleed area but vary the radial 
position of the effective bleed area (i.e., if 
bleed was desired at an inner radius, the outer 
holes were covered, or vice versa.), 

since the bleed enun was fitted with pressure 
taps, flow velocities could be attained by 
securing the suction to the plenum and recording 
the pressure drop in the plenum due to the flow 
over the bleed holes. By covering all but one 
row of holes the velocity at a desired radius 
could be obtained by using the isentropic relation 
for pressure and Mach number or Bernoulli's 


equation. 


F. EXPERIMENTAL PROCEDURE 


Initially six nozzle rings were tested in the apparatus 


illustrated in figure (7). Stagnation pressure applied to 


the nozzle stack, controlled by means of a regulator valve, 


was increased to a maximum of 130 pounds per square inch. 


Boundary layer bleed was attained by connecting the vacuum 


pump directly to the bleed plenums. The flow was observed 


using shadowgraph and Schlieren techniques. A typical run 


consisted of establishing the sensitivity of the Schlieren 
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D 


(or shadowgraph) system to the desired level, starting the 
vacuum pump and allowing it to reach its lowest pressure, 
and then increasing the pressure to designated levels where 
pnotographs were taken. 

The number of nozzles was then changed from six to four 
and the spacing of the nozzle rings (the throat height) was 
increased from .05 inches to .09 inches. This change 
decreased the exit Mach number from 4.5 to 3.9. Since the 
total pressure loss across the shock was less, the back 
pressure did not need to be lowered to below atmospheric. 
Therefore На ejector nozzles were removed, Removing the 
ejectors allowed for longer runs at high pressures due to 
a decrease of approximately forty percent in the mass flow 
from the reservoir. 

The four nozzle rings were tested in the apparatus 
modified as shown in figure (32) in a fashion similar to that 
previously described for six nozzles. One importent differ- 
ence was the use of a 200 cubic foot tank to increase the 
bleed capability of the vacuum pump. The pump was connected 
to the tank and was able to lower the pressure inside it to 
approximately one-half pound per square inch. The bleed 
plenums were then connected to the bleed tank and the suction 
was controlled with a gate valve. 

Flow conditions were observed and photographed at 
pressures up to the maximum attainable pressure of 250 pounds 


per square inch. Photographs were taken with and without 
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boundary layer bleed to compare the resulting flow fields. 
The bleed hole pattern was also altered to obtain pressures 
in the bleed plenum due to the velocity of the flow across 
the holes at various radii as discussed in Section III-E, 

A pointed probe was inserted into the flow at various 
locations and different total pressures to determine regions 


of subsonic and supersonic flow. 
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ШЕН ГИС ФАН) CONCLUSIONS 


А, RESULTS 

As was previously mentioned, six nozzle rings were 
initially installed in the test apparatus shown in figure 
(32). Normal shocks were obtained in all except the end 
nozzles, but attempts at moving the shocks beyond the nozzle 
exit to form a "barrel" shock were unsuccessful at pressures 
of up to 130 psi." Figures (33)(a) and (b) are Schlieren 
photographs of the six nozzles at pressures of 70 and 130 
psi, respectively. Before attempting higher pressures it 
was decided to increase the throat area in order to lower 
the exit Mach number and thus make it easier to move the 
shocks out of the nozzles. Since the sum of the throat areas 
must remain smaller than the feed area to prevent choking at 
the ends of the stack, the number of nozzle rings was reduced 
to four and the throat height increased to .09 inches. This 
resulted in an exit Mach number of 3.9 and increased the Mach 
number at the feed holes to .9. 

Three distinct flow patterns were observed in the flow 
over the range of pressures up to a maximum of 250 psi. 
Figures (34), (35) and (38) are Schlieren photographs 
showing the development of the first stage flow structure. 


Figure (37) is a schematic of the flow just prior to its 


* # 
Applied pressures are gege pressures unless otherwise 
specified. 
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(a) 


(b) 130 PSI, 


FIGURE (33) SCHLIEREN PHOTOGRAPHS OF NORMAL SHOCKS 
STANDING IN SIX NOZZLES 
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(a) NO FLOW 


FIGURE (34) SCHLIEREN PHOTOGRAPHS OF FOUR NOZZLES 
BEFORE "STARTING" 
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FIGURE (35) SCHLIEREN PHOTOGRAPHS OF CENTER THREE 
NOZZLES STARTING 
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(a) Flow at Pressure Less Than 50 PSI. 
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(b) Resulting Sonda ler attached tó Nozzle Wall 


FIGURE (36) SCHEMATIC OF COANDA JET FORMATION 
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SCHLIEREN PHOTOGRAPH OF FLOW JUST BEFORE 


FIGURE (38) 


FLOW STRUCTURE 


TRANSITION TO SECOND 
TOTAL PRESSURE 125 PSI. 
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transition to the second stage at 145 psi. It can be 
readily observed that failure of the end nozzles (actually 
"half" nozzles) to start complicates the flow structure a 
great deal. Coanda jets appear to be attached to end nozzle 
walls creating a "bubble" of relatively slow flow. A slip 
layer results between this slow region and the Coanda jet 
and appears in the Schlieren photographs as a light area. 
It was originally considered that these bright regions could 
be expansion waves turning the flow back to each wall. 
Inserting a probe into this region as shown in figure (42-b) 
established that the flow was in fact subsonic since it 
failed to produce oblique shock waves from the tip of the 
probe. Pressure readings were also recorded on the wall at 
the nozzle exit and by assuming isentropic flow from there 
to ambient pressures (calculating upstream), the Mach number 
was calculated and found to be subsonic. Table (5) lists 
the pressure data taken and the calculated Mach numbers. 
Figures (39) and (40) show the pressure and Mach number at 
the wall as functions of total pressure. The slip-line 
causes a bright area in the Schlieren photographs due to the 
higher entropy and temperature generated from friction, 
leading to an area of lower density. 

Careful analysis of the photographs of the flow up to 
pressures of approximately 145 psi. (e.g. before transition 
to the second phase) show a shock wave starting normal to 


the wall within the end nozzles and then curving to the 
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Table 5 Wall Static Pressure at Nozzle Exit and 
Calculated Mach Number. 


Total Gage Wall Static” Mach Number” 
Pressure (psi) Pressure (psi) at Wall 
0 12,94 „43 
20 12.42 „50 
40 Js 259 
60 IELO . 66 
80 10.51 20 
100 deer 79 
120 9.28 . 84 
130 - / 8.79 - 789 
135 -- 7 12207 - /. „90 
140 82.56. 7.10.63 о 
145 LOCO ШО З Hoo 59 
150 10.69 / 10:78 20907 2268 
155 10,79 68 
160 10278 .68 


ж 
The second values listed resulted from decreasing the total 
pressure from 150 psi (illustrating the hysteresis effect). 
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angle of the nozzle wall, where the slip layer forms 

between the Coanda jet and the subsonic "bubble." The slip 
layer continues to curve back toward the end wall. Observable 
in the lower jet are two small, nearly normal shock waves — 
one where the shock attached to the wall begins to curve 
outward and one from the trailing edge of the nozzle ring. 
These combine to slow the flow to subsonic velocities as 
testing with the probe revealed. The jet on the upper 
nozzle is smaller and faster which is attributable to a 
slightly smaller opening at the throat. Figure (38) shows 
the alternately light and dark regions of a series of shock 
diamonds in the upper Coanda jet, indicating a higher Mach 
number than in the lower nozzle. The slip line in the upper 
case also bends back toward the wall. 

The high velocity jets occur for one of two reasons. 
First, since the end nozzles are really only "half-nozzles," 
the throats are only half the thickness of those in the 
center, in this case on the order of .045 inches. With 
these small openings located at the ends of the nozzle stack 
where the flow is entering axially at near sonic conditions, 
it has a difficult time "turning the corner" to exit out the 
first throat opening. If this were the case, whatever flow 
did make the 90° turn may attach itself to the nozzle wall 


as shown below, 





fir separation 
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A second, more feasible explanation is that regardless of 
whether or not the flow has trouble turning the corner, the 
flow that travels through the throat starts as a jet of air 
near the wall at low pressures (in the order of 50 psi.) is 
shown in figure (36-a). This jet entrains the air naar the 
nozzle wall and thus lowers the pressure in the wall region. 
The jet of air then moves toward the wall and attaches to 
it, displaying the "Coanda effect." As the flow in the jet 
moves Out radially it slows and increases in pressure, 
forcing the slip line back toward the wall and creating a 
"bubble," observable in figure (35). As the total pressure 
increases toward the transition point of 145 psi., the 
bubble becomes longer as the flow in the jet takes a longer 
distance to reach a higher pressure. 

The three nozzles in the center established normal 
shocks quite readily and moved outward to the nozzle exit 
with increasing pressure. As can be seen in the Schlieren 
photograph of figure (38), connecting these three waves 
would result in a cylindrical shock bowed inward at the 
center, re > states that a loss of total pressure 
in an axial direction will manifest itself in a shock wave 
that is conical, tapering inward in the direction of 
decreasing total pressure, rather than strictly cylindrical. 
For a stack of nozzle rings fed from both ends, as Was the 
case in this experiment, such a loss in total pressure would 
result in the condition observed at pressures between 90 and 


ШІ БӨЗ. chat as, a tapering inward, from both ends. Total 
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pressures different from the region mentioned, however, did 
not result in a distortion of this type (see figures (41) 
through (48)). Seemingly, if axial total pressure losses 
were observable in a concave shock at pressures near 100 psi. 
they would also be observable at much higher pressures 
rather than disappear near 145 psi. 

At a pressure of approximately ten atmospheres, transition 
to the second phase flow structure occurred. Figures (41) (a) 
and (b) are Schlieren photographs of the flow at total 
pressures of 150 and 175 psi. This transition consisted of 
a Sudden rise in the wall pressure, and therefore a drop in 
the Mach number, as shown in figures (39) and (40), respec- 
tively. The Coanda jets no longer form a subsonic bubble, 
but instead the slip layer extends upward at an angle of 
19 degrees. The cause of the transition to the condition 
shown in figure (43) is the movement of the normal shock 
segments beyond the trailing edges of the nozzles. Since 
the flow leaves the nozzles at approximately a five to ten 
degree angle, obligue waves bend the ends of the normal 
shock segments outward giving them a concave shape. Shock 
diamonds form behind the oblique waves serving to return 
the flow to subsonic conditions.  AVCO Everett in Ref. 2 
predicts that disturbances traveling downstream from the 
trailing edges would only affect the cylindrical shock wave 
locally, causing outward protrusions on the surface of the 


shock. Though the oblique shocks are rather large 
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(a) 150 PSI. 





MAS PSI. 





FIGURE (41) SCHLIEREN PHOTOGRAPHS SHOWING CURVED SHOCK 
SEGMENTS BEYOND TRAILING EDGES OF NOZZLES 
AFTER TRANSITION TO SECOND PHASE 
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(a) 150 PSI. - SUBSONIC 
BEHIND CURVED SHOCK 
SEGMENT 


HI 150 PSI. - SUBSONIC 
IN WALL REGION 


FIGURE (42) SCHLIEREN PHOTOGRAPHS WITH PROBE SHOWING 
REGIONS OF SUBSONIC FLOW 





В СЗ 


SHOCK 
Se ae 


SLIP-LINE 
ISIS JET = — 
7 


М <1 


К 


NS SOS 
w 

WEAK 

Es SHOCK 


N DIAMONDS 
= 


ES 





SLIP-LINE 


EL DIAMONDS 
М<1 


И oe d E 
KIN Е IN 


FIGURE (43) SCHEMATIC OF THE FLOW STRUCTURE ATA 
TOTAL PRESSURE OF 175 PSI, 
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FIGURE (44)  SCHLIEREN PHOTOGRAPHS OF FLOW JUST BEFORE 
DEVELOPING INTO THIRD STAGE — TOTAL 
PRESSURE OF 200 PST. 
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(a) 220 PSI. 





"240 PSI. 





FIGURE (45) SCHLIEREN PHOTOGRAPHS SHOWING DEVELOPMENT 
OF FLOW TO THIRD STAGE 


3, 





ee 


"disturbances," this prediction appears to hold by the fact 
that the ends of the normal shock segments "protrude" 
outward. Flow exiting in a more nearly horizontal direction 
would minimize the strength of these oblique waves, thereby 
making the shock segments less curved and closer to joining 
into a complete cylindrical shock wave. 

The second phase is also characterized by the disappear- 
ance of the subsonic "bubble" in each end nozzle, As the 
normal shock in the passage above (or below) the end nozzle 
moves beyond the trailing edge, a pressure mismatch occurs 
between the very low pressure region behind the shock and 
the Coanda jet coming off the trailing edge on the opposite 
side of the nozzle ring. This causes the Coanda jet to be 
diverted upward (or downward), and the resulting slip-line 
to extend upward (or downward). The angle that the slip-line 
forms increases from 19 degrees at transition (150 psi.) to 
29 degrees at the maximum pressure of 250 psi. This change 
is due to the increase in the pressure mismatch created at 
the trailing edge. Since the trailing edge is blunt, the 
Coanda effect will continue to force the jet to higher angles 
as it tries to move around the trailing edge. This is not 
the case with a sharp trailing edge where the angle of the 
Coanda jet would remain the same, 

AS the pressure was increased from 150 psi. to 200 psi. 
the flow maintained the basic configuration shown in figure 


(43). Figures (41), (42) and (44) are Schlieren photographs 
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of the flow in this range Of pressures. At 200 psi the 
oblique waves from the trailing edges are longer, better 
defined, and at a smaller angle. The curved shock segments 
have moved to approximately a quarter of an inch beyond the 
trailing edges. At approximately 210 psi. the flow structure 
is altered a second time. 

Transition to the third stage, unlike the abrupt change 
at 145 psi.,was a gradual development of normal shock segments 
in the nozzle wakes. Between pressures of 205 and 215 psi. 
the normal segments slowly appear, along with two weak oblique 
shocks at the trailing edges where one strong wave had been 
at lower pressures. Increasing the total pressure to 250 psi. 
changed the flow structure only by moving both the normal and 
curved segments of the shock slightly farther out from the 
trailing edges. Figure (47) is a schematic of the fully 
developed flow. Figure (46) is a Schlieren photograph of 
the flow at the maximum pressure applied, and shows that the 
normal portions of the shock are at a distance of approximately 
three eighths to a half of an inch. This distance is short 
of creating a lasing cavity wide enough to allow the theoret- 
ically available power in the flow to be extracted, which for 
four nozzles is 1.7 inches. For a spacing of the nozzles 
в, (s)/(r,-r,) т ОЕ ете т. is the shock location. 
As a fraction of the minimum shock location allowable to 
extract the power available, (r Lo), the location of the 


m 


shock wave is: 
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FIGURE 


(46) 





SCHLIEREN PHOTOGRAPH OF FLOW AT MAXIMUM 
TOTAL PRESSURE OF 250 PSI. NORMAL SHOCK 
SEGMENTS APPROXIMATELY 3/8 INCH FROM 
TRAILING EDGES 
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FIGURE (47) SCHEMATIC OF THE FLOW STRUCTURE AT A 
TOTACTPRESSURE OF 250 PSI. 
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Boundary layer bleed was applied over the range of 
pressures with little if any effect on the flow. Figures 
(48) (a) and (b) are Schlieren photographs of the flow at 
200 and 250 psi, respectively with a one-half psi. bleed 
pressure applied to the first three rows ОҒ bleed holes on 
the upper plate. At a total pressure of 250 psi the flow 
entrained between the slip-lines is diverted toward the 
bleed plate. This was the only noticeable change caused 
by the bleed over the entire range of pressures. Attempts 
were made at applying enough suction to start the end 
nozzles but were unsuccessful. 

Homogeneity of the flow in the lasing cavity is a 
necessity for an operational laser. This was not the case 
for the fully developed flow at 250 psi. since the end 
nozzles failed to accelerate the flow to supersonic condi- 
tions, Other disturbances such as the oblique shock waves 
from the trailing edges should not be present to produce 


homogeneous conditions in the lasing cavity. 
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(a) 200 PSI. 


(b) 250 PSI. 


FIGURE (48) SCHLIEREN PHOTOGRAPHS SHOWING FLOW WITH 
BOUNDARY LAYER BLEED APPLIED TO UPPER PLATE 
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В. CONCLUSIONS AND RECOMMENDATIONS 
Becsrenoso:stheseylindrical nozzle stack at a total 
pressure of 250 psi. established a shock wave approximately 

a half of an inch from the nozzle exit. The shock wave 
consisted of two segments normal to the flow (actually 
shock "bands" around the nozzle stack) connected by shock 
segments curved inward. Failure of the end "half-nozzles" 
to start resulted in subsonic flow near the end walls for 
the entire range of total pressures investigated. This 
prevented the normal shock from reaching the end walls. 
Since subsonic flow existed at the end walls and the 
shock did not reach the walls, the use of boundary layer 
bleed had very little visible effect on the flow structure. 

Flow visualization became very distorted at high 

pressures due to the extreme density gradients and the 
cylindrical geometry of the shock waves. 

The following recommendations are made for further 

research on this subject: 

1. Increase the diameter of the nozzle rings to allow 
for more mass flow to enter the nozzle stack without 
decreasing structural integrity and without choking 
in the feed passages. 

2. Decrease the thickness of the nozzle rings to lower 
the expansion ratio and thus lower the exit Mach 
number, 


3. Ensure sharp trailing edges on the nozzle rings, 
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To help avoid starting problems in the end nozzles, 
use half of a nozzle ring (sliced along the edge) 

at each end. This would eliminate the "half-nozzle" 
arrangement and the Coanda jets. 

Investigate total pressures in excess of 250 psi. 
Install pressure taps on the center nozzle wall to 
determine Mach number. Optical distortion due to 
the large density gradients and the cylindrical 
geometry results in difficulty observing shock 


angles from a probe inserted in the flow. 
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APPENDIX A 


THEORETICAL LASING POWER AVAILABLE 


Intensity is defined as radiation power per unit area, 
with units of kilowatts per square foot. Assuming a 
Gaussian intensity distribution on each mirror as shown 


in figure (A-1) the intensity is 
_ 2 

I(r) = I e i 
To satisfy the boundary condition of 

E EE в/о) 
Е = Хх, then 

Que DAY : 
uneretrore: 

I(r) = І e : (A-1) 
Equation (A-1) yields a Gaussian peak centered at r=0. To 


obtain a ring with Gaussian cross section having the peak 


intensity at r-r, define 
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FIGURE (A-1) LASING MIRROR INTENSITY LOADING 





Inserting this into equation (A-1) yields 





ыз с 


I(r) + Т е E "Y kwatts/cm* 2 


To obtain the total power output, the intensity must be 


integrated over the mirror area. 








т 2 
m = 
Power = f f I. exp(=t V2 1% r ar de 
= P г -г 
е 0 е 
Em 
Deren]? ar (A-2) 
Р г EEE 


An alternate method of finding the power available from the 


flow is 


(specific) 
energy 


Power (mass flow rate) 


Power = E (kjoules/1b,) m (1b,/Sec) 
P= Em kwatts 


or 


P = dp, a Lë 3 28 : (5257 
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Equating equations (A-2) and (A-3): 





r-r 2 T 
I y r expl --7/2 ] dr - rS b0gag (05/05) Lë Bel MP р 
(А-4) 
where the only unknown in equation (A-4) is the outer radius 
OF Che mirror, En 


In solving the integral Q on the left side of equation 


(A-4) define 








r “Y 
making 
\/2 
ах = ах 
r Y 
and 
ro-rY 
r= лү 





Making the corresponding changes in the limits of 


integration: 
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the integral now takes the form 


-V2 к-ті Е 
г ch o; 














Or 
Q » FCR ак ак + жете; T < ax 
= E e x x e 
gd v2 yz (А-5) 


Since the second term on the right side of equation (A-5) 
contains an integral of an odd function between symmetric 


limits, the integral goes to zero leaving 


NA ONE Ur eT ахи. (A-6) 


Multiplying and dividing the right side of equation (A-6) by 


"T /2 and rearranging, 


"/2т r (r-r ) v2 2 2 
Q = = J — == Өте ах 
4 o ут 


(A-7) 


Noting that the integral in equation (A-7) is the definition 


of the error function gives 


О = \т/2 Y (r-r.) erf (V2 ) : (A-8) 
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Substituting equation (A-8) for the integral in equation 


(А-4): 
г 101/2 Е (For) erf(V2 )] = кра (р/р,) (а/а )М Е , 


SE, Solving for Y , 


(./2) * 1/2? z^ *C (A-9) 


Y = 
where 
az 4 Ге L Mo E 
^/т/2 т erf (y2 ) 
Since 


А] 
1 
Y 


г = r + C > (A-10) 


The values in the expression for C are: 


ro = 1.5 inches 
ро = (.002378) (32.174) (10) = .765 16 /ЕЕ? 
ааа 1119 £t/sec 
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4 


> 


(РИРо) м=4.5 = 01745 


Е = 5 kwatts/ (lb, /sec) 


3 


I = 3 Б = 2.8 x10 kwatts/£t” 


р 
егҒ(/2 ) = 2(.9773-.500) = .9546 . 


Using these values in the expression for C: 
С = ,2692(L) inches” . 


For DM 1.5 inches, equation (A-10) becomes 


г - 7/2251 2692" L 
m 


where L is the length of the nozzle stack in inches. 
Dimension L is illustrated in figure (7). 

As an example the power available from six nozzles was 
calculated. The theoretical lasing power available was 


found to be 
Im 
r-r 


J | г expl VE 1% ar 
Por r -r 


e e 
271, er£(y2 ) dau? [(F(F-r_)] 


21(3) (30.52 /144) (.9546) (п/2) ?[1.635(.135)] 





И 
N 
= 
H 


Power 


1 


И 


Power = 32.16 kilowatts. 


| 
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АРРЕМОТХ В 


CALCULATIONS FOR DETERMINING INLET (FEED) AREA 


Figure (13) shows the dimensions of the nozzle ring. 


The area of the hole segment between two centers is 


A = length x width 


= =» (hole diameter) 


= Em (1/8) x (3/8) 
2 


А = „1407 in К 


Since there are four segments: 


I 


Az 4(.1407) 


2 


A 26627 1n 777 


1 


The area of the remaining half circles is 


БД 1,m(3/8)° 


psi 


ны 
| 


A, — quje in* y 


Therefore the total inlet area is 


ы 
lI 


in A1 Ad 


E m 
Bin = 1.104 in 
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APPENDIX C 


DETAILS OF THROAT HEIGHT CALCULATIONS 


The feed flow and the flow exiting through the nozzles 


must satifsy the continuity equation: 


И 
пм е 


nca "nozzle 


== XITXAX = 
:م2‎ іп іп 6p*U*A (С-1) 
where the factor 2 in equation (C-1) results from the 
nozzle stack being fed from both ends, 
Using the equation of state, and noting that velocity 
is the product of Mach number and the speed of sound, 


equation (C-1) can be written 








Solving for the throat height: 


P: 
Е L in E 
ЕЕ К э. c 
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: ; L 
Multiplying the right side by (5-2 апа (ж—) * , 
O O 
equation (C-2) can be written as 
P T 
Hl no Ol 0,1% | 
h 7 (—) М. (-5-) (ву) (—) (—) e (C23) 
бт ın Ро P To Tig 


Using the isentropic equation relating stagnation temperature 
and pressure to static temperature and presure respectively, 
and applying them to the Mach equals one condition, 
equation (C-3) becomes 

T 


за 1 in 1 O 
h = (==) М m (у) (93353) лит: 


1 
2 


Or 


м. (P. /P.) 
в = ‚0917 م ك‎ (С-4) 
В 


Noting that 


a | 
A и (ү-1) „2. Y/Y-1 


equation (C-4) reduces to 


" 1 
l = s. 0917 Min Е 


(T/Tin) ° (To/Tin 
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or 


Ytl 
A ee 
0917 м. (ТИТ)? | | 


(ү-1) 
.0917 eege м, 


209 7 М. 


м 


(1+.2 Ше ) 
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APPENDIX D 


` AXIAL VARIATION OF MACH NUMBER AND PRESSURE 


From Ref. 4, influence coefficients for flow in a 


channel are given as 





| ү-1 „2 
am? — 1+ М) an 
zu. (D-1) 
м2 1-м2 
апа 
ар м^ ад 
T ii e Di (D-2) 
1-м 


assuming isentropic flow and constant specific heat and 
molecular weight. Using the model shown in figure (D-1) 


and an inlet area of 1.105 square inches, 





Bout E Il 

А. 16 
т 

АО б (ey (Aa) 
out 16 in 

E 2 52 

А = .06906 іп 2 

out 
Since 
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NEN 06906 
Ay 1.105 - „06906 
or 
QA oa 
KG ES 06667 o 


Substituting this value into equations (D-1) and (D-2) and 


rearranging yields 


ү-1 2 
M.. =M. -M í RE (.06667)] (D-3) 
itl i 3 2 S 
1 - M. 
i 
and 
үм,“ ‚ 
i 


A program was written for the Hewlett Packard model 9830-A 
computer to calculate the axial variation of Mach number and 
pressure using equations (D-3) and (D-4). The program 
determines Mach MU EE and pressure at stations one through 
thirteen using values of inlet Mach number and total 
pressure (PSF) as inputs. A listing of the program AXIAL 


and the results of the program are found in this appendix. 





Table D-1 Symbols for Program AXIAL 


Computer Notation Algebraic 
M1 M; A 
PØ ры 


ШЕ) 


PÈ ) 
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Description 


Inlet Mach number 


Stagnation 
pressure (PSF) 


Mach number at 
subscripted 
station 


Static pressure 
at subscripted 
station 





COMPUTER PROGRAM "AXIAL" 


LISTING OF 
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RESULTS OF PROGRAM 
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APPENDIX E 


DETERMINING POINTS ON THE NOZZLE WALL 
——————— 1 - 


Following a development similar to Ref. 4, but changing 
the characteristic network as shown in figure (E-1), 
equations can be developed to determine points along a 
right running characteristic proceeding in the - & direction. 


Along the characteristic segment 


3» 3 ` sin 6 

Iu oan S сіп S= dn (E-1) 
on E 

2 2 

and 

3 3 | 

y 221000 СП л ети iW ПО qz (Е-2) 

1 95 1 5 


where the negative sign on the right side of equation (E-2) 
results from traveling in the negative -& direction when 
proceeding from point 1 to point 3, Assuming a very small 
characteristic mesh, the integrands on the righe sides of 
equations (E-1) and (E-2) can be taken to be constant, 

With this assumption and some rearranging, the equations 


become 


біп 0 біп 6 
Sn нә K о 





ER 
+ 5 [-sin ЕТ ЕА 
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CHARACTERISTIC 


FIGURE (E=1) 








LEFT-RUNNING 
CHARACTERISTIC 


STREAMLINE 


MODIFIED CHARACTERISTIC ELEMENT 
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1 sin 0 
+ > [-sin На Е, 


AE13 - sin Ы, 
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AEPENDIX E 


А. EJECTOR NOZZLE DESIGN 

To assure starting of the cyiindrical nozzles a back 
pressure must exist which is lower than the total pressure 
downstream of the shock, For a stagnation pressure of ten 
atmospheres, the total pressure downstream of a normal 
shock at a Mach number of five is less than one atmosphere. 
Ejectors were designed to lower the back pressure to a value 
below the total pressure downstream of the shock assuming 
Mach-five flow. 

The models used for finding the flow conditions at the 
ejector nozzles are shown in figures (F-1) and (F-2). 


Applying the continuity equation 
DL EE аи © E 


the momentum equation 


2 2 = 2 Т 
mu ро в OU AF PA; y (5-2) 


and the energy equation 


2 2 n 2 
Az (GP, +U,%/2) + р ОА (СТО +0, (оу = p,U3A, (GT, +U, ED) 
(F-3) 
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FIGURE (F-1) EJECTOR NOZZLE MODEL 
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BECURE ел) ТӘРГУТЕМ OF EJECTOR FLOW MODEL 
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to the control volume shown will yield the flow conditions 
necessary at station one to assure the desired back 
| pressure. 
The following assumptions are made: 
(1) isentropic flow from station 4 to 2 


(2) isentropic flow from reservoir to ejector nozzle exit 


(3) Pz = 1 atm. 
(4) Pry > GCG E. 
(5) Р, = Р, 


(6) Тт, = 520°R 
(7) M = 5.0 at shock 


Using normal shock relations the known flow at station (5) 
can be used to obtain conditions at station (4). For 


М = 5,0, the pressure 


O 


P, = (29) (39.99) 
Pow „548 atm. , 
density 
04/5 = 5,0 


ра 7 (5.0) (o3) (05/09) 
- (5.0)(.02378 g) (.01134) 


Е поради ее | 
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and temperature 


T¿/ To —— 5G 


H 
ll 


so тт 


(Gores) (520) (21667) 


| 


СИТ СУЕ 


FH 
И 


can easily be found. 


Using the isentropic relations 
P ZP) = (p,./p ) = (T,/T)) |, 


all conditions at station two can be found in terms of the 


static pressure: 


Po = na (PAI 


aa S 
2 


ll 


P2 
and 


аа ае 


2 


H 
Il 


Since ті. = My ВЕБЕ velocity at 2 can also be found: 


m, = Pog (2A5) = mA 


Ee? 





where the area is taken as (2:A,) since the control volume 
inlet area (A^) is only one-half of the total exit area of 


the flow.  Rearranging: 


m, 
U, = рә А» ЕА. х (144) 
(6.463) (144) 
Po T o 
55,9 
U, = =a : 
2 Po 


All conditions at station two are now known in terms of Ро. 
Again using isentropic relations, conditions at station one 


can be determined in terms of Ру: 


Р 
| T -(y-1)/Y 
а: 
1 T4 Py 
-.4/1,4 
> в" : 
T4 >20 (5) ор 
апа E 
ра + Pp, (——) 
1 0 Ра 
-1/1.4 2 
"EN EUN >; lb/ft? 
1 Р, m 
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To find the velocity at station one: 






u vet 2 
Tr, T, (1 +—— M; DEL 72 
2 u, 
Since Му - YRTY , 
T4 2YRT4 
U, = Le: DIG g ft/sec . 


Hence all conditions at station one are known in terms of 
Pi. 
bie Р, (=P,) is known, equations (E-1) and (E-2) can be 
combined to find U, and 05. A computer program was written 
to vary the values of Dy and р, and to iteratively increase 
Р, from a set value until it satisfied the energy equation. 
An explanation of the program EJECTOR, a listing, and the 
results in tabular form are found in this appendix. Figure 
(26) shows the results in graphical form. The broken line 
at a pressure of 0.617 represents the total pressure down- 
stream of the shock assuming Ро (upstream) = 10 atm. and 


М = 5.0. Acceptable combinations of ру and Dz must result 


in a pressure P, below this line. 
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Bee EXHAUST LENGTH 

The length of the exhaust stack must be long enough to 
ensure complete mixing of the primary flow from the ejector 
nozzles and the secondary flow from the cylindrical nozzles. 

The model used for mixing of the two flows is shown 
merigure (E-3). 

Bailey [Ref. 6] defines a non-dimensional spread 


parameter which is a good approximation to experimental data: 
U 
с = 24 2 (F-4) 


where 


Unean Е 2 i 


The non-dimensional similarity variable for the mixing 


zone was defined to be 

Dis o E (F-5) 
A value of n = 1.84 gives the outer edge of the mixing zone 
where the velocity equals 99% of the centerline value. 


The minimum value of the mixing length is where the edge 


of the mixing zone touches the side walls, or when 
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FIGURE (F-3)  EJECTOR FLOW MIXING MODEL 
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Substituting Y = and equation (F-4) into equation (F-5), 


the minimum mixing length, Хы сап Бе found: 





35914 
X m Сї ах 
min n 
Ber en (D4 = ру) 
Ui 2(1,84) 
(U, +U.) 
= 6.522(4.85 - 1.75) 
2 
л 
„522 
ү 


Using values from program EJECTOR for the combination of 


Di and D4 chosen: 


1860 ft/sec 


Ui 


Во 


664 ft/sec 


Putting these values in equation (F-6): 


004922 664 
Eu Л тш > 


13.72 inches .. 


>, ae 
min 


The exhaust stack must therefore extend a minimum of 13,72 
inches beyond the ejector nozzle exit. Figure (F-4) shows 
the length of the exhaust stack with the ejector nozzles 
installed. 
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GC. NOZZLE SHAPE 


For the chosen combination of Dy = 1.75 inches and 


Dz = 6.5 inches 


Uy = 1860 ft/sec . 


(See section D of this Appendix.) This velocity is clearly 


supersonic. To determine the throat area required for the 


given exit diameter and flow velocity, the exit Mach number 


most first be found. Since 


жор Us 
m 

L u“ 
yRT : 


and 


= E ee 
T = ЗЕ ои Je в 


then, after substituting and rearranging, the Mach number 


can be expressed as 


where 
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Inserting the appropriate values: 


For isentropic flow through the nozzle 


ү+1 
1 2 Y-l „2 2 (y-T) 
re] 


A/A* 


For M = 2.49 


aa 


Selving for throat height: 





a A M 
4 2.61 
m PAS 
А* = ,921 in ; 
or 
4A* 
х = 
р ТІ 
D* = 1,08 inches, 


The ejector nozzles were fabricated from plexiglas with the 
throat diameter of 1.08 inches. Figure (F-4) shows the 


details of the ejector nozzle geometry. 
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FIGURE (F-4) EJECTOR NOZZLE GEOMETRY 








Ша EJECTOR 
The program EJECTOR solves the continuity and momentum 
equations for different combinations of nozzle diameter and 
exhaust stack width assuming a starting value of pressure 
at stations one and two of 0.37 atmospheres. The pressure 
Pi is increased in increments of 0.0001 atmospheres a 
the flow conditions at 1, 2 and 3 satisfy the energy equation. 
When the energy equation is satisfied, the values are printed, 
and the process is repeated with new values of D, and Da. 
Inputs to the program include conditions immediately 
downstream of the shock: 0, В Е : Ta (ORUM. 
Pa (atm.) and m, ua pa | ‚ and total temperature (°R) 
and total pressure (atm.) at station one. 
Outputs are D, and D4 (both in inches), Py ib) in 


atmospheres, and Ui and Us, Беси in ft/sec: 
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Computer Notation 


Table F-1 


А А2 АЗ 


Р1 
Rl 


Ml 


P2 P3 P4 
R2 R3 R4 


M2 M3 M4 


B2 T3 T4 


02 U3 


Algebraic 


А) А А5 


DI 


A 


Баби и ve 
Ша Mo ща m, 
ы 
R 


e mme тт. 


Symbols for Program EJECTOR 


Description 

Flow area (in^) 

Static pressure (atm.) 
Static density (UPS EE) 
Mass flow (1b /sec) 
Specific heat ratio 


Gas constant for air 
(59 3 ft-lbe/lb, ор) 


Specific heat € constant 
pressure 


776:26 ft-lb./BTU 

Static temperatures (?R) 
Stagnation temperature (?R) 
Stagnation pressure (atm.) 


Velocities at stations 1, 2 
and 3 respectively (ft/sec) 


Diameter of One of the Ewo 
ejector nozzles (inches) 


Length of one side of exhaust 
stack. 











D. LIST OF COMPUTER PROGRAM "EJECTOR" AND RESULTS 
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